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Abstract 

We present an effective model for the generic behaviour of hadron masses 
and phase shifts at finite temperature which shares basic features with recent 
developments within the PNJL model for correlations in quark matter. On 
this basis we obtain the transition between a hadron resonance gas phase 
and the quark gluon plasma in the spirit of the generalized Beth-Uhlcnbeck 
approach where the Mott dissociation of hadrons is encoded in the hadronic 
phase shifts. Here we restrict ourselves to low-lying hadronic channels and 
perform a discussion of recent lattice QCD thermodynamics results from 
this perspective. We find agreement in the asymptotic regions while for the 
description of the transition itself the inclusion of further hadronic channels as 
well as a selfconsistent determination of the continuum thresholds is required. 


1 Introduction 

Simulations of lattice QCD (LQCD) are in practice the only reliable ap¬ 
proach to QCD thermodynamics which covers the broad region of strongly 
interacting matter properties from the hadron gas at low temperatures to 
a deconfined quark gluon plasma phase high temperatures. Recently, finite 
temperature LQCD simulations have overcome the difficulties of reaching the 
low physical light quark masses and approaching the continuum limit which 
makes this theoretical laboratory now a benchmark for modelling QCD under 
extreme conditions [1.2). In spite of this, effective models as such prove valu¬ 
able for interpreting LQCD results as well as for extensions of QCD based 
ab-initio knowledge to the regions of the QCD phase diagram which are inac¬ 
cessible to LQCD and high energy heavy ion collision experiments. These are 
the regions high net baryon number and isospin densities at very low temper¬ 
atures as probed by observations of astrophysical objects like compact stars 
as well as in neutron star mergers and supernovae collisions. 

The aim of this paper is to construct a combined effective model re¬ 
producing the two known asymptotics of the finite temperature equation of 
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state: the hadron resonance gas at low temperatures and the free quark-glnon 
plasma with perturbative corrections in accordance with the recent LQCD 
simulations. This endeavour started about ten years ago with an ansatz 
for temperature- and mass dependent width in a hadron spectral function 
|3> 4J. As was then demonstrated in J5] the nice description of LQCD data 
at that time stemmed from the fact that low-lying meson states below a 
mass threshold of ~ 1 GeV were not subject to spectral broadening so that 
at high-temperatures where the majority of states of the hadron resonance 
gas disappeared nevertheless a Stefan-Boltzmann limit was attained which 
accidentally coincided with the LQCD one because the number of low-lying 
hadronic degrees of freedom matches that of quarks and gluons. Note that 
the idea to use a spectral function approach to describe the dissociation of 
the hadrons in a resonance gas with increasing temperature has recently also 
been formulated in Ref. [B], but due to a lacking dynamical model of hadrons 
as bound states of quarks and gluons, the thermodynamics of the latter had 
there to be added by construction from a comparison with LQCD data. 

The next step [5] was to lift the lower mass threshold and let all hadrons 
be subject to spectral broadening leading to a decrease of their contribution 
to the total pressure and to add the pressure of the quark-gluon system 
according to the meanheld treatment of the Polyakov-loop generalized NJL 
model [7J. Some of the arbitrariness of the choice of the temperature- and 
mass dependence of the spectral width T(T, M ) was removed in [8j where the 
latter was taken from a recently developed approach to chemical freeze-out 
[9] • This model could also make very nice predictions for the parton fractions 
in the temperature region of the crossover transition from hadronic to quark- 
gluon matter, in accordance with, e.g., Ref. [10]. However, this model did 
still not yet provide a microscopic description of the temperature- and state- 
dependent spectral broadening of hadrons in the medium. 

In the present contribution we want to join these ideas for the description 
of a Mott-hadron resonance gas (MHRG) with a Beth-Uhlenbeck approach 
m to pion and sigma meson dissociation ra in its recent generalization 
from the NJL model H3 to the PNJL model case jT4l IT5l ITB] . The latter 
approach is based on the in-medium modification of (quark-antiquark-) scat¬ 
tering phase shifts which encode the spectral properties of the complex-valued 
meson propagators. For an alternative but equivalent formulation, see mi. 
In order to join the Beth-Uhlenbeck formulation of Mott dissociation with 
the hadron resonance gas, we will postulate a generic behaviour of the phase 
shifts in selected hadronic channels which are temperature dependent and 
embody the main consequence of chiral symmetry restoration in the quark 
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sector: the lowering of the thresholds for the two- and three-qnark scatter¬ 
ing state continuous spectrum which triggers the transformation of hadronic 
bound states to resonances in the scattering continuum. The phase shift 
model is in accordance with the Levinson theorem (see, e.g., Refs. [181 [19]) 
which results in the vanishing of hadronic contributions to the thermodynam¬ 
ics at high temperatures m- While in this work we will restrict ourselves to 
a few states only, the model is of course applicable to the whole spectrum of 
hadronic states in the particle data book, which is presently under way [20]. 

2 Quark and hadron mass spectrum 

The underlying quark and gluon thermodynamics is described within a Nf = 
2 + 1 PNJL model at the mean-field level 

Ppnjl(T) = P cond (T) + P FG (T) + U(T), (1) 

with the Nf = 2 + 1 parameters determined analogously to [210 resulting 
in a dynamically generated vacuum light quark mass of m( 0) = 394.7 MeV, 
and a strange quark mass of m s (0) = 604.745 MeV for the parameters A 
= 600 MeV, G s A 2 = 2.403, m 0 = 5.4285 MeV and m s = 136.459 MeV. 
One obtains the temperature dependence of both the light quark mass m(T ) 
and the strange quark mass m s (T) as a solution of the corresponding gap 
equations. 

We make the simplifying ansatz that the hadron masses are constant and 
temperature independent up to their Mott temperature, where they hit the 
two-quark threshold m t h r ,M(P) (three-quark threshold m t h r ,B(7 1 )) for mesons 
(baryons), respectively. After that temperature, we assume a linear rise of 
the resonance mass with temperature, being identified with the resonance 
width T i (P) , unique for all hadrons 

M r (T) = Mi( 0) + R(T), Tj(T) — (T — T Mott f)6(T - (2) 

where for the constant slope we choose a = 2.5. The Mott temperatures 
TMott.i can be determined as soon as the temperature dependence of the light 
quark mass m{T) is known. They are determined from the condition 


M i (T Mot t,i) 

?Whr,i(pMott,i)> 

(3) 

"hhr,M(P) 

= (:2 - N s )m(T) + N s m s (T ), 

(4) 

rrhhr,B(T) 

= (3 - N s )m(T) + N s m s (T ), 

(5) 


2 The values in |21j are to be corrected, see http: //3fcs .pendicular .net/psolver 
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Figure 1: Temperature dependent hadron mass spectrum used in the model. 


where N s = 0,1, 2 for mesons {i = M ) and Ns = 0,..., 3 for baryons {i = B). 

The mass spectrum ([3]) determines the corresponding Mott dissociation tem¬ 
peratures as TMott,i[MeV] = 303, 338, 196, 252 (for hadrons i = 7r,K,p,N, 
resp.) after which the hadrons became unbound state, see Fig. [lj 

3 Mott effect and temperature dependent hadron 
phase shifts 

We define temperature dependent phase shifts 6i(s;T ) for each hadronic 
channel (species) i as functions of the Mandelstam variable s = M 2 by a 
generic behaviour which holds for both types of hadrons, mesons and baryons 
and encodes the following characteristic features 

• a bound state with mass M; is recognized by a jump of the phase shift 
by +7T at s = Mf ; 

• above the continuum threshold s t h r ,i = mj( hri the phase shift drops to 
zero which is reached at s maXii = s t h r ,i + NfN 2 where N t = 2 for mesons 
and Ni = 3 for baryons; 
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• For T > TMott.i the step function of the true bound state gets smeared 
out over a width Fj(T) resembling a resonance in the continuum. 

Our ansatz is given by 

Si(s-,T) = 

+ 

where the auxiliary function F(s) = sin(s/r 2 )@(r 2 7r/2 — s) + 0(s — r 2 7r/2) 
has been introduced in order to ensure that the phase shift at s = 0 shall 
always be zero, even at higher temperatures, where large values of the width 
parameter in the Breit-Wigner like ansatz would otherwise spoil this con¬ 
straint. A microscopic calculation of the finite temperature pion phase shift 
which shows the correct behaviour has been performed in [13]. It has been 
used to constrain the above ansatz for F(s). 


, ( s-M?(T) \ 

F(s) —I- arctan - - —— —— 

- m 2 hr i )0(m 2 hr4 + iV 2 A 2 - s) 


hhr,i 


m thr,i + ^A 2 ~ S 

N? A 2 


b (6) 



Figure 2: Generic meson phase shift as defined by this model, here for the 
pion (N t — 2,i — 7r) at six temperatures between 150 and 400 MeV. 

The behaviour of the phase shift ([6]) as a function of the squared mass 
variable s is shown for the example of the pion in Fig. [2] The chosen values 
of six temperatures in steps of 50 MeV, T = 150... (50)... 400 MeV cover 
the region of the Mott transition temperature T Mo tt, 7 r where the pion bound 
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state (jump of S n (s) from zero to tt at s = M 2 ) merges with the edge of the 
continuum of scattering states at s — m t \ r)r . At this temperature the mass 
gap (binding energy) vanishes which in Fig. [2] is recognizeable as the plateau 
at S n = tt and consequently for higher temperatures T > Tm 0 u,tt the pion 
phase shift starts at zero without reaching the value tt again. For higher 
temperatures the flattening of the hadron phase shift describes the vanishing 
of the hadronic correlation from the thermodynamics of the MHRG system. 


4 Mott hadron resonance gas (MHRG) 


For the thermodynamic potential, i.e. the pressure, of two-particle correla¬ 
tions in quark matter we have derived |T2] the formula 


P,(T) 


di 

di 


d 3 p [°° doj 


(2vr) 3 J 

30 dpp 2 


TT 

' ds 


1 


27r 2 Jo 27r a Jp 2 + s 


+ s)Si(s-,T), 


(7) 


which we apply here to mesons and baryons, with di being the degeneracy of 
the state i and in the distribution function = l/[exp(o;/T) =F 1], where 
the upper (lower) sign holds for mesons (baryons). One can use s = M 2 and 
integrate over the mass variable M instead of s, and changing the order of 
integrations to obtain 




dpp 2 r°° dM M 


27 r 2 Jo tt a Jp 1 + M 2 

By partial integration over M we obtain 


fi (Vp 2 + AT%(M 2 ;T). (8) 


PAT, = Td, f°° ^ r dM rind + e-vTVTP/T) 
Jo 2tt- Jo 


1 d5i(M 2 ;T) 


TT 


dM 


(9) 


the generalized Beth-Uhlcnbeck formula for the partial pressure of species 
% in the MHRG. When the phase shifts degenerate to a step function, i.e., 
when m{T) —> oo, T^T) —» 0, so that 5j(M 2 ; T) = tt6(M — Mj), we have 


1 d5i(M 2 \T) 
tt dM 


6{M-Mi), 


so that the M-integration becomes trivial and gives 


( 10 ) 


P,(T) = T d. 




(11) 
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the relativistic ideal quantum gas for a particle species i with mass Mj. For 
the total pressure of the model, we have 

P(T) = Y, Pi( T ) + ^pnjl(T) + P P eit(T) , (12) 

i=M,B 

where the pressure of the quarks and gluons ([Tj) , and a perturbative pressure 
correction have been added. The latter results from 0(a s ) interactions of 
quarks and gluons in a plasma for momenta exceeding the cutoff A to which 
the nonperturbative interactions in the PNJL model are restricted. Details 
are described in Refs. M- Note that at variance with these references, we 
use here a temperature dependent, regularized running coupling [22, [23] 


r 2 — c 2 

where r = 3.2T, c = 350 MeV and N c = Nf = 3. 



a(T) = 


12vr 


1UV C - 2Nf \ln(r 2 /c 2 ) 




Figure 3: Temperature dependence of the pressure for linear scale (left panel) 
and logarithmic scale (right panel) |Tj. 

In Fig. [3] we compare the total pressure of the MHRG model (black solid 
line) with LQCD data from Ref. [T] (star symbols) and show also the partial 
pressure contributions from all components in the linear (left panel) and 
logarithmic scale (right panel). With the present schematic model, we can 
quantify both asymptotic limits, the pion gas limit at low temperatures T < 
100 MeV and the quark-gluon plasma limit at high temperatures T > 300 
MeV. In the region between these temperatures we observe a discrepancy 
between the model and LQCD data which is attributed to mainly two missing 
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aspects: i) the inclusion of higher lying hadronic resonances and ii) the back 
reaction of the MHRG on the temperature dependence of the quark masses 
and thus on the continuum thresholds which trigger the onset of the Mott 
dissociation of hadrons. 


5 Conclusions 

In this work an effective model is constructed which is capable of reproducing 
basic physical characteristics of the hadron resonance gas at low temperatures 
and embody the crucial effect of hadron dissociation by the Mott effect. The 
generalized Beth-Uhlenbeck form of the partial pressures is constructed for 
each hadronic channel. Numerical results show that the simplifying ansatz 
for the temperature dependence of both, the mass spectrum and the phase 
shifts of hadronic channels give results in qualitative agreement with recent 
ones from LQCD [TJ. The differences at the present level of sophistication of 
our effective model are well understood and basically two conclusions for its 
further development can be drawn. 

First, one can include higher lying hadronic states into the model. Second, 
and most importantly, the back reaction of the hadron resonance gas on the 
chiral condensate and thus to the quark masses and the continuum thresholds 
which derive from it must be taken into account. This latter point will lower 
the onset of chiral restoration and hadron dissociation towards the result from 
lattice QCD. It shall further pronounce the character of a smooth crossover 
transition between MHRG and QGP by quantifying the relative contributions 
of hadronic and partonic degrees of freedom (partial fractions) to the QCD 
thermodynamics. 
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